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Abstract: In the context of climate change, precipitation is predicted to become more intense at the 
global scale. Such change may alter soil microbial communities and the microbially mediated carbon and 
nitrogen dynamics. In this study, we experimentally repackaged precipitation patterns during the growing 
season (from June to September) of 2012 in a semi-arid temperate steppe of the Xilin River Basin in 
Inner Mongolia of China, based on the 60-year growing season precipitation data. Specifically, we 
manipulated a total amount of 240 mm precipitation to experimental plots by taking the following 
treatments: (1) P6 (6 extreme precipitation events, near the 1* percentile); (2) P10 (10 extreme 
precipitation events, near the 5" percentile); (3) P16 (16 moderate precipitation events, near the 50% 
percentile); and (4) P24 (24 events, 60-year average precipitation, near the 50" percentile). At the end of 
the growing season, we analyzed soil microbial community structure and biomass, bacterial abundance, 
fungal abundance and bacterial composition, by using phospholipid fatty acid (PLFA), real-time 
quantitative polymerase chain reaction (RT-qPCR) and 16S rRNA gene clone library methods. The 
extreme precipitation events did not change soil microbial community structure (represented by the ratio 
of PLFA concentration in fungi to PLFA concentration in bacteria, and the ratio of PLFA concentration 
in gram-positive bacterial biomass to PLFA concentration in gram-negative bacterial biomass). However, 
the extreme precipitation events significantly increased soil microbial activity (represented by soil 
microbial biomass nitrogen and soil bacterial 16S rRNA gene copy numbers). Soil fungal community 
showed no significant response to precipitation events. According to the redundancy analysis, both soil 
microbial biomass nitrogen and soil ammonium nitrogen (NH4-N) were found to be significant in shaping 
soil microbial community. Acidobacteria, Actinobacteria and Proteobacteria were the dominant phyla in 
soil bacterial composition, and responded differently to the extreme precipitation events. Based on the 
results, we concluded that the extreme precipitation events altered the overall soil microbial activity, but 
did not impact how the processes would occur, since soil microbial community structure remained 
unchanged. 
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1 Introduction 


Shifts in precipitation regime can be predicted by general atmospheric circulation models (GCMs), 
with the expectation of increasing extreme conditions, such as precipitation events with large 
amount but low frequency (Easterling et al., 2000; Weltzin et al., 2003; Pachauri et al., 2014). 
This change will be a remarkable contrast to the current precipitation events with small amount 
but high frequency distributed across relatively short intervals that have been historically defined 
in the semi-arid ecosystems (Heisler-White et al., 2008; Hu et al., 2017). Many studies have 
concentrated on the influences of precipitation pulse, rewetting, drought and shifting amounts of 
precipitation on ecosystems (Landesman and Dighton, 2010; Hueso et al., 2012; Bouskill et al., 
2013; Sorensen et al., 2013). Nevertheless, little is investigated on the effects of shifts in 
precipitation events on ecosystems (in particular on an ecosystem like a semi-arid steppe), 
especially on the belowground environment, including soil microbial community biomass and 
structure (Knapp et al., 2008; Sorensen et al., 2013; Knapp et al., 2017). 

Persistent shifts in precipitation patterns that fall into the range of statistical extreme climate 
events can impose substantial effects on soil nutrient cycling and soil microbial community of the 
grassland ecosystems (Fierer and Schimel, 2002; Knapp et al., 2002). Soil microbial community 
structure and biomass are sensitive to shifts in precipitation patterns (Steenwerth et al., 2005; 
Sorensen et al., 2013; Canarini et al., 2016). However, soil bacteria and fungi response differently 
to precipitation patterns (Barnard et al., 2013). Soil fungi are thought to be more resistant to water 
variability than soil bacteria due to fungal abundant mycorrhiza through which water can be 
transferred from one point to another (Allen, 2011; Bell et al., 2014; Zhang et al., 2016). Soil 
bacteria are more sensitive to drought and rewetting at both DNA-based and RNA-based levels in 
California grasslands (Barnard et al., 2013) due to their cell structure and life strategy (Che et al., 
2015; Kaisermann et al., 2015). A previous study suggested that in California grasslands, it is the 
extreme desiccation and rewetting event significantly affects the present and potential soil active 
bacterial community but not the fungal community (Barnard et al., 2013). Alternatively, in 
laboratory experiments, increasing the number of rewetting events not only decreases soil CO2 
loss but also changes soil microbial biomass dynamics (Fierer and Schimel, 2002). 
Drying-rewetting cycling experiment can decrease soil bacterial biomass while increase the ratio 
of fungal to bacterial growth (Bapiri et al., 2010). There are some other studies still suggest that 
seasonal precipitation variability alters soil fungal community abundance and consequently 
influences soil nutrient cycling (Cregger et al., 2012; Bell et al., 2014). Shifts in frequency of 
precipitation events may also affect the aeration status of soil matrix and the soil functional 
microbes (directly and indirectly), and influence ecosystem processes more than gradual climate 
change does (Horz et al., 2004, 2005; Jentsch et al., 2007; Zhang et al., 2013). Despite these 
generalities, in semi-arid ecosystems, extreme precipitation events do not impose on soil nutrient 
availability and subsequently on soil microbial community in the same way as drying-rewetting 
events or precipitation amount will do (Belnap et al., 2004; Heisler-White et al., 2008; Chen et al., 
2015). 

In this study, we established a field experiment by simulating different precipitation patterns 
(frequency and amount) in a semi-arid steppe in the Xilin River Basin, China. Our overall goal 
was to ascertain whether shifts in precipitation patterns, predicted by the GCMs, would cause 
changes in the structure and biomass of soil microbial community and subsequently affect soil 
physical-chemical properties such as soil pH, soil carbon and soil nitrogen. We hypothesized that 
the extreme precipitation events (large in amount but low in frequency) would decrease the 
biomass of soil bacteria and fungi and change the structure of soil microbial community. 
Moreover, we further addressed the following questions in this study: would the extreme 
precipitation events affect the biomass and abundance of soil bacteria and fungi as well as the 
structure of soil microbial community? If so, would soil microbial community structure correlate 
with soil nutrient pools? 
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2 Materials and methods 


2.1 Study area 


The field experiment was conducted at the Inner Mongolia Grassland Ecosystem Research Station, 
Chinese Academy of Sciences (43°32'N, 116°40’E; 1200 m a.s.1), located in a semi-arid steppe in 
the Xilin River Basin, Inner Mongolia, China. The research field has been fenced off since 1979. 
The region is characterized by a temperate continental climate with an annual mean temperature 
of —0.48°C and a mean annual precipitation of 358 mm (Hao et al., 2010). Precipitation during the 
growing season highly accounts for 89% of the total annual precipitation, of which 75% (about 
240 mm) is considered as the ecologically effective precipitation (EP; recorded daily 
precipitation >3 mm) (Hao et al., 2011, 2012, 2013). Soil cation exchange capacity is 15.7 
cmol/kg and soil pH value is 7.5. Vegetation community in this steppe comprises a mixture of 
annual grasses, annual and perennial forbs, and perennial shrubs. The xeric rhizomatous grass 
(Leymus chinensis) and needlegrass (Stipa grandis) are the dominant species (Liu et al., 2017). 
The soil is loamy sand or dark chestnut. 


2.2 Experimental design 


Precipitation treatments were carried out during the growing season (from 1 June to 30 September) 
of 2012. The primary aim of these treatments was to alter the dry interval duration and the 
frequency of precipitation events. It should be noted that we maintained the constant precipitation 
of 240 mm during the experimental period, approximately equivalent to the 60-year average 
during the growing season. Four precipitation treatments were done in 16 plots in total through a 
random block design, with four replicates for each treatment. The four treatments included: (1) P6 
(6 extreme precipitation events, near the 1* percentile); (2) P10 (10 extreme precipitation events, 
near the 5" percentile); (3) P16 (16 moderate precipitation events, near the 50" percentile); and (4) 
P24 (24 events, 60-year average precipitation, near the 50™ percentile) (Fig. 1). Precipitation 
intervals were 20-, 12-, 7- and 5-d for the four treatments, respectively. The amount of each 
precipitation event was not constant and was determined according to the 60-year average over 
the same season. Average precipitation amount was 10.0 (ranging from 3.8 to 21.8 mm), 15.0 
(ranging from 6.0 to 26.5 mm), 24.0 (ranging from 11.6 to 41.5 mm) and 40.0 mm (ranging from 
22.2 to 72.6 mm) for P24, P16, P10 and P6 treatments during the experimental period, 
respectively. The sprinkling can was used to simulate the precipitation scenario. Based on a 
long-term data analyzed by Huang et al. (2010), we applied larger events at 2—3-d intervals to 
guarantee that the maximum daily precipitation in the experimental plots would never take more 
than 24.0 mm. Volumetric soil water content (SWC) was monitored every week using a TDR 300 
Soil Moisture Meter (Spectrum Technologies, USA) at depths of 5 and 20 cm. 

Each plot was settled as an area of 2 mx2 m, and was surrounded by a metal flashing extending 
approximately 10 cm above and 40 cm below the ground to isolate roots and to hold back the 
lateral water flux. A rain-exclusion shelter (3 mx3 m), consisting of a steel frame supporting a 
clear 0.8-mm thick fiberglass reinforced polyester roof with only minor shading effects (90% 
light transmission), was installed on 1 June in each plot to eliminate ambient precipitation. The 
shelters covered the plots during the experimental period. To evaluate microclimate effects, we 
measured the photosynthetic active radiation (PAR) using a LI-190SB quantum sensor (LI-COR, 
Inc., Lincoln, NE, USA), and the air temperature using a HMP45C temperature probe (VAISALA, 
Woburn, MA, USA) under the shelters with insuring that they were measured in an open space 
close to the plots. We found no significant differences in PAR and air temperature between plots 
and ambient surroundings. 

The longest available historical precipitation data during the period 1953—2010 were collected 
from the Meteorological Administration of Xilin Gol League. However, only the EP data were 
used in our data analyses. 
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Fig. 1 Occurrence probability of the number of precipitation events during the growing season of 2012 based on 
an estimated probability function calculated from the 60-year average precipitation of the growing season. Four 
precipitation treatments were designed: (1) P6 (6 extreme precipitation events, near the 1* percentile); (2) P10 (10 
extreme precipitation events, near the 5" percentile); (3) P16 (16 moderate precipitation events, near the 50" 
percentile); and (4) P24 (24 events, 60-year average precipitation, near the 50" percentile). 


2.3 Soil sampling and analysis 


At the end of the experimental period (30 September), three soil cores (diameter of 3 cm and 
depth of 10 cm) were taken in the first three replicate plots for each precipitation treatment, and 
then they were mixed after elimination of plant litter and roots. Soil samples were sieved through 
a 2-mm mesh to avoid diluting the influence of surface litter on soil microbial community 
composition. Then, each soil sample was divided into two parts. One part was maintained fresh 
for measurements of the composition and biomass of soil microbial community, and soil 
inorganic nitrogen. The other part was air-dried to determine soil organic carbon (SOC), soil pH 
and soil total nitrogen (TN). 

SOC was determined using the H2SO4-K2Cr207 oxidation method (Islam and Weil, 1998). TN 
was analyzed by the Kjeldahl digestion method (Mckenzie and Wallace, 1954). Inorganic nitrogen 
containing ammonium nitrogen (NH4-N) and nitrate nitrogen (NO3-N) was extracted with 2 M 
KCI and analyzed by a continuous-flow ion auto-analyzer (AutoAnalyzer 3, Seal Analytical, 
England). Soil pH was measured in a 1:2.5 soil:water suspension. 

Soil microbial biomass carbon and microbial biomass nitrogen were measured using the 
fumigation-extraction method (Brookes et al., 1985; Vance et al., 1987) with several 
modifications. In brief, 10 g of paired soil samples were weighed and then incubated at 4°C for 24 
h. One sample was fumigated for 24 h with ethanol free CHCl; away from light, while another 
sample was not fumigated. At the end of the fumigation, both fumigated and non-fumigated soil 
samples were extracted from 0.5 M K2SOg in a shaker for 30 min. The extractions were passed 
through the filter paper after 30 min standing and were determined using a vario TOC analyzer 
(Elementar, Germany). 


2.4 Phospholipid fatty acid (PLFA) analysis 


In this study, PLFA was used to evaluate coarse-scale changes in soil microbial community 
structure following the study of Frostegard et al. (1993). Briefly, 8 g of sampled fresh soil at the 
end of experiment was extracted with a single-phase mixture of citrate buffer, methanol and 
chloroform on a horizontal shaker at room temperature for 2 h. Phospholipids were quantified by 
a gas chromatography-mass spectrometry (GC-MS, Agilent 6850, USA). The fatty acids (14:0, 
al4:0, al5:0, 115:0, 116:0, 16:l@5c, 16:1@7c, i117:0, cy17:0, a17:0, 18:1@8, 18:1@7c, cy19:0 and 
cy20:0) were chosen to present the bacterial biomass (Frostegard et al., 1993). The 18:26, 9 was 
taken to indicate the fungal biomass (Xue et al., 2005). The 10Me18:0, 10Me16:0 and 10Me17:0 
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were used to determine the Actinobacteria. The branched, saturated fatty acids (114:0, al5:0, 115:0, 
116:0, i17:0 and al7:0) were used as a classification of gram-positive bacterial biomass (Zogg et 
al., 1997). The cyclopropane-unsaturated and monoenoic fatty acids (16:l@5c, 16:1@7c, 16:1@9c, 
17:l@8c, 18:l@5c, cy17:0 and cy19:0) were used as a classification of gram-negative bacterial 
biomass (Xue et al., 2005). All results were converted to a dry weight basis. 


2.5 Real-time quantitative polymerase chain reaction (RT-qPCR) 


DNA extracted from the sampled fresh soil (0.35 g) at the end of experiment was isolated using 
the soil DNA extraction kit (Mobio Powersoil DNA Isolation Kit, Carlsbad, USA). Then, it was 
estimated with a Nano Drop ND-1000 UV-Vis Spectrophotometer (Nano Drop Technologies, 
USA) following its manual. 

The RT-qPCR was performed with specific primers on the StepOne Plus Real-time PCR System 
(Applied Biosystems, USA). The RT-qPCR reactions were performed in 20 uL volume mixture 
containing 10 uL SYBR® Premix Ex Taq™ (Takara Bio Inc., Japan), 10 uM each primer and 2 uL 
template DNA solution. Forward primer 338F (5'-ACT CCT ACG GGA GGC AG-3') and reverse 
primer 518R (5'-ATT ACC GCG GCT GCT GG-3') were used to identify the population size of 
total bacteria (Edwards et al., 1989). Forward primer gITS7 (5'-GTG ART CAT CGA RTC TTT 
G-3') and reverse primer ITS4 (5'-TCCTCCGCTTATTGATATGC-3') were used to amplify for 
fungal ITS region. The RI-qPCR reactions used to analyze bacterial 16S rRNA gene were 
performed as following conditions: denaturation 5 min at 94°C, followed by 42 cycles of 94°C (60 
s), 56°C (45 s) and 72°C (45 s), and then final extension 10 min at 72°C; and the reactions used to 
analyze ITS gene were performed as follows: denaturation 5 min at 95°C, followed by 40 cycles of 
95°C (20 s), 54°C (25 s) and 72°C (45 s), and then final extension at 72°C for 10 min (Ihrmark et 
al., 2012). Finally, melting curve analyses were performed to ensure the specificity of the PCR 
(polymerase chain reaction) amplicon (Bustin et al., 2009; Zhang et al., 2017). 

Standard curves were produced by a serial dilution of sequences extracted from a plasmid 
(pGEM-T easy vector, Promega) containing 16S rRNA gene and ITS gene fragments amplified 
from Escherichia coli DH5a competent cells (Takara Bio Inc., Japan). They were linear 
correlation over six orders of magnitude (R?>0.98). 


2.6 16S rRNA gene clone library 


In order to analyze the soil bacterial community composition in surface layer (0—10 cm depth), we 
constructed a bacterial 16S rRNA gene clone library for the four treatments in this study. 
Amplification of approximately 500-bp 16S rDNA was performed in a MyCycler thermocycler 
(Bio-Rad, Hercules, CA, USA) with forward primer 27F (5'-GAGTTTGATCMTGGCTCAG-3') 
and reverse primer 519R (5'-TATTACCGCCGCKGCTG-3') (Hill et al., 2003). The PCR reactions 
were performed as following conditions: 5 min for initial denaturation at 94°C, 35 cycles of 94°C 
(45 s), 56°C (45 s) and 72°C (50 s), and then final extension 10 min at 72°C. For the PCR 
products, gel electrophoresis analyses were performed in 1% (wt/vol) agarose gel and then the 
products were recovered using a gel purification kit (Axygen, USA). The obtained PCR products 
were ligated into a pUCm-T vector (Promega, USA) and transformed into Escherichia coli DH5a 
competent cells. The positive clones were screened on indicator plates by color-based clone 
selection. Furthermore, the positive clones were confirmed with PCR reactions using specific 
primers T7 (5'-TAATACGACTCACTATAGGG-3') and M13 (5'-CAGGAAACAGCTATGACC-3’) 
of pUCm-T vector, and the PCR reactions were performed the same as 16S rRNA gene 
amplification. 


2.7 Data analysis 


The data of PLFA, microbial abundance and soil properties were analyzed using one-way 
ANOVA that identified treatments as fixed effect in SAS 8.0 (Institute Inc., Cary, USA). 
Significant levels for all statistical tests were set at P<0.05 level. Before the significant tests, the 
normality of error terms was calculated by the Kolmogorov-Smirnov test for goodness of fit, and 
the homoscedasticity was evaluated using the Levene test for equality of variances. NMDS 
(Non-Metric Multi-Dimensional Scaling) for PLFA of soil microbe was also performed to 
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determine how soil microbial community structure was related to precipitation events. Stress 
<0.05 in NMDS was taken for representing a good reflection of the overall community structure. 
RDA (redundancy analysis) was performed to test the relationship between PLFA sample and 
environmental variables. NMDS and RDA were performed by the R 3.2.5 (R Development Core 
Team, 2016) with vegan package. 


3 Results 


3.1 Variations of volumetric soil water content (SWC) 


Variations of volumetric SWC under different precipitation treatments during the growing season 
of 2012 are shown in Figure 2. Generally, the precipitation treatments altered SWC significantly 
at the 20 cm depth (F=3.32, P=0.02), but the influence was not significant at the 5 cm depth 
(F=2.27, P=0.08). The SWC under the four precipitation treatments showed significant 
differences on the same day (Figs. 2a and c). The mean SWC was higher under the P24 treatment 
than under the other three treatments during the whole experimental period at the two soil depths 
(Figs. 2b and d). The mean value of SWC varied from 11.5% under the P6 treatment to 14.4% 
under the P24 treatment at the 5 cm soil depth (Fig. 2b), and from 14.4% under the P6 treatment 
to 16.1% under the P24 treatment at the 20 cm soil depth (Fig. 2d). 
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Fig. 2 Variations of volumetric soil water content (SWC) at 5 (a, b) and 20 cm (c, d) soil depths under four 
precipitation treatments during the growing season of 2012. * in the left panel means significant differences of 
SWC among the treatments on the same day at P<0.05 level. Different lowercase letters in the right panel indicate 
significant differences of seasonal average SWC among the treatments at P<0.05 level. Bars mean standard 
errors. 
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3.2 Soil properties and microbial biomass 


During the growing season of 2012, the precipitation treatments exhibited no significant influence 
on SOC (F=2.80, P=0.13), TN (F=0.00, P=0.96) and soil pH (F=0.29, P=0.60) in all plots (Table 
1). However, these treatments significantly affected the instantaneous soil properties, i.e., NH4-N 
(F=7.77, P=0.03) and NO3-N (F=6.94, P=0.03). There was a large range of NH4-N values 
(4.36-10.43 mg/kg) among the four treatments, with the values higher in the extreme 
precipitation events (P6 and P10 treatments) than in the other precipitation events (P16 and P24 
treatments; Table 1). The values of NO3-N ranged from 0.88 to 1.54 mg/kg across the four 
treatments. The lowest NO3-N value was found under the P6 treatment, which was significantly 
lower than that under the P24 treatment (Table 1). 


Table 1 Effects of precipitation treatments on soil physical-chemical properties 


Treatment SOC (%) TN (%) NH,-N (mg/kg) NO3-N (mg/kg) pH 
P6 2.09+0.07 0.26+0.01 7.774£1.51" 0.88+0.14° 7.3140.12 
P10 2.07+0.01 0.25+0.01 10.43+0.27° 1.3540.22” 7.3340.28 
P16 1.80+0.07 0.24+0.01 7.2441.25° 0.91+0.12° 7.40+0.05 
P24 2.05+0.07 0.25+0.01 4.36+0.89° 1.54+0.19* 7.55+0.09 


Note: Four precipitation treatments were designed: (1) P6 (6 extreme precipitation events, near the 1* percentile); (2) P10 (10 extreme 
precipitation events, near the 5" percentile); (3) P16 (16 moderate precipitation events, near the 50" percentile); and (4) P24 (24 events, 
60-year average precipitation, near the 50" percentile). SOC, soil organic carbon; TN, total nitrogen; NH4-N, ammonium nitrogen; 
NO;-N, nitrate nitrogen. Different lowercase letters indicate significant differences among treatments at P<0.05 level. Mean+SE; n=3. 


The extreme precipitation events of the P10 treatment significantly influenced soil microbial 
biomass nitrogen (F=18.89, P<0.01; Fig. 3a). Soil microbial biomass nitrogen under the P10 
treatment (167.51 (+6.61) mg/kg) was twice higher than that under the P24 treatment (71.69 
(£0.51) mg/kg). However, the extreme precipitation events (P6 and P10 treatments) did not 
significantly alter soil microbial biomass carbon (F=0.66, P=0.44; Fig. 3b) and the ratio of 
microbial biomass carbon to microbial biomass nitrogen (F=0.39, P=0.55; Fig. 3c). 
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Fig. 3 Variations of microbial biomass nitrogen (a), microbial biomass carbon (b), and the ratio of microbial 
biomass carbon to microbial biomass nitrogen (c) under different precipitation treatments. MBN, microbial 
biomass nitrogen; MBC, microbial biomass carbon. Different lowercase letters indicate significant differences 
among the treatments at P<0.05 level. Bars mean standard errors. 


3.3 Response of soil microbial community structure to precipitation events 


Response of soil microbial community structure to different precipitation treatments is shown in 
Figure 4. There were no significant influences of precipitation treatments on the PLFA 
concentration of aerobic bacteria (F=0.37, P=0.56; Fig. 4a), anaerobic bacteria (F=0.03, P=0.86; 
Fig. 4b), Actinobacteria (F=0.00, P=0.96; Fig. 4c), fungi (F=1.71, P=0.23; Fig. 4d) and total 
bacteria (F=0.53, P=0.49; Fig. 4e). However, all of the PLFA concentrations were lower under the 
P6 treatment than under the P24 treatment, with an exception of Actinobacteria. The lower values 
of PLFA concentration in aerobic bacteria and anaerobic bacteria were found under the P6 and 
P10 treatments. 
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Fig. 4 Variations of phospholipid fatty acid (PLFA) concentration of aerobic bacteria (a), anaerobic bacteria (b), 
Actinobacteria (c), fangi (d), and total bacteria (e), and variations of the ratio of PLFA concentration in fungi to 
PLFA concentration in bacteria (Fungi/Bacteria; f) and the ratio of PLFA concentration in gram-negative bacterial 
biomass to PLFA concentration in gram-positive bacterial biomass (G-/G+; g) under different precipitation 
treatments. Bars mean standard errors (n=3). 


The extreme precipitation events did not significantly affect the ratio of PLFA concentration in 
fungi to PLFA concentration in bacteria (F=1.70, P=0.24; Fig. 4f), with the ratio values ranging 
from 0.11 to 0.13. The maximum value occurred under the P24 treatment (0.12). In addition, the 
extreme precipitation events also did not significantly affect the ratio of PLFA concentration in 
gram-negative bacterial biomass to PLFA concentration in gram-positive bacterial biomass 
(F=0.24, P=0.63; Fig. 4g). However, as shown in the NMDS plot (Fig. 5), soil microbial 
community structure characterized by the overall PLFA concentration differed between the 
extreme and mean precipitation events (stress=0.015). 


3.4 Relationship between soil nutrient and microbial community structure 


We found that soil NH4-N (F=3.8, P=0.05) and microbial biomass nitrogen (F=5.4, P=0.03) were 
significant in shaping soil microbial community structure in RDA. The first and second RDA axes 
explained 53.3% and 1.4% of the total variation in soil microbial community structure, 
respectively (Fig. 6). 

3.5 Response of soil bacterial and fungal abundance to precipitation events 


The changed precipitation patterns significantly influenced soil bacterial 16S rRNA gene copy 
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Fig. 5 Non-Metric Multi-Dimensional Scaling (NMDS) ordinations of soil microbial community structure 
assessed by PLFA analysis under the P6, P10, P16 and P24 precipitation treatments (stress=0.015) 
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Fig.6 Biplot of PLFA sample-environment variables using the redundancy analysis (RDA) 


numbers (F=9.05, P=0.02; Fig. 7a). The extreme precipitation events tended to increase soil 
bacterial abundance. Soil bacterial 16S rRNA gene copy numbers varied from 9.38x10° to 
13.20x10° copies/g dry soil under the four precipitation treatments (Fig. 7a), with the maximum 
value under the P10 treatment and the minimum value under the P24 treatment. However, 
precipitation treatments showed no significant impact on soil fungal abundance (F=0.60, P=0.63; 
Fig. 7b). Soil fungal ITS gene copy numbers ranged from 4.19x10° to 4.9610? copies/g dry soil 
under the four precipitation treatments, with the maximum value under the P6 treatment and the 
minimum value under the P10 treatment. Soil fungal ITS gene copy numbers were much lower 
than soil bacterial 16S rRNA gene copy numbers (Figs. 7a and b). 


3.6 Response of soil bacterial composition to precipitation events 


The identified 16S rDNA sequences of soil bacteria contained abundant major bacterial phyla 
under the four different precipitation treatments. The dominant phyla were Acidobacteria, 
Actinobacteria and Proteobacteria (Fig. 7c). Compared with the P24 treatment, the mean relative 
abundance of Actinobacteria increased by 21.0% and 43.8% under the P10 and P6 treatments, 
respectively, while the mean relative abundance of Acidobacteria decreased by 49.4% and 67.4% 
under the P10 and P6 treatments, respectively. We also found other common phyla, such as 
Gemmatimonadetes, Planctomycetes, Verrucomicrobia and Planctomycetes, under the four 
treatments. However, Firmicutes was found only under the P6 treatment. 
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under different precipitation treatments, and stacked bar charts of variations in soil bacterial composition at the 
phylum level under different precipitation treatments (c). Different lowercase letters indicate significant 
differences among the treatments at P<0.05 level. Bars mean standard errors (n=3). 


4 Discussion 


The extreme precipitation events decreased the frequency of precipitation but increased the 
intensity (amount) of precipitation, therefore extending drought intervals (Gong et al., 2004; 
Heisler-White et al., 2008; Pachauri et al., 2014). These events are different from the patterns of 
drought, dry-rewetting and precipitation amount (Berard et al., 2011; Barnard et al., 2013; Zhang 
et al., 2013). The extreme precipitation events included the integrated dynamics of three 
simulating experiments: drought, dry-rewetting and precipitation amount. We found that extreme 
precipitation events significantly decreased SWC, especially at the 5 cm soil depth. Such change 
may alter soil microbial community and the microbially mediated carbon and nitrogen dynamics. 
Soil microbial community in the temperate steppe of the Xilin River Basin showed a higher 
variability responding to extreme precipitation events during the growing season. We 
hypothesized that an integrated effect of dry intervals, magnitude of precipitation and more 
dry-rewetting events would modify soil microbial community structure toward a higher ratio of 
fungi to bacteria and a lower bacterial abundance. This pattern would be related to the availability 
of soil nutrient, as well as SWC (Placella et al., 2012; Bell et al., 2014; Gschwendtner et al., 
2014). However, we found that extreme precipitation events did not alter soil fungal abundance, 
the ratio of soil microbial biomass carbon to soil microbial biomass nitrogen and the ratio of 
PLFA concentration in fungi to PLFA concentration in bacteria, but increased soil bacterial 
abundance. Furthermore, the higher NH4-N content may have been responsible for the higher 
bacterial abundance in the extreme precipitation events. These findings indicated that the extreme 
precipitation events (P6 and P10 treatments) did not modify soil microbial properties during the 
growing season. Meanwhile, soil fungal community showed a high resistance to extreme 
precipitation events. In contrast, the extreme precipitation events affected soil bacterial abundance 
and significantly influence soil nutrient availability, especially for soil inorganic nitrogen. Similar 
results have been found by Bell et al. (2014) in Chihuahuan Desert grasslands that 25% changes 
in seasonal precipitation frequency could significantly affect soil nutrient and microbial 
properties. 

Changes in the amount, frequency and timing of precipitation may impact soil carbon and 
nutrient pools (Nielsen and Ball, 2015; Wilcox et al., 2015). The extreme precipitation events did 
not change SOC and TN, but altered soil nutrient availability including NH4-N and NO3-N 
contents, and then increased soil microbial activity (represented by soil microbial biomass 
nitrogen). What's more, the changed NHg4-N and soil microbial biomass nitrogen were 
significantly correlated with soil microbial community in the extreme precipitation events. This 
correlation was likely linked to an alteration of soil microbial heterotrophic respiration and 
rhizosphere activity (Chen et al., 2009; Li et al., 2012), which feedback to the nutrient balance 
and cycling of underground ecosystem (Broeckling et al., 2008). The changes in soil inorganic 
nitrogen reflected more instantaneous responses of soil nutrient bank in the extreme precipitation 
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events. Results from previous studies had also revealed that precipitation pulse stimulated 
mineralization and nitrification of soil nitrogen (Dijkstra et al., 2012), while drought influenced 
nitrogen losses through denitrification (Evans and Burke, 2013). 

Fungi are generally considered more resistant to extreme climate than bacteria, especially under 
drought conditions (Gordon et al., 2008; de Vries et al., 2012; Li et al., 2012). Similarly, we also 
found that soil fungal abundance and biomass kept stable both at gene and community levels. The 
changes in the ratio of microbial biomass carbon to microbial biomass nitrogen and the ratio of 
PLFA concentration in fungi to PLFA concentration in bacteria imply a transformation in soil 
microbial community structure (Zeglin et al., 2013). However, we did not find that the ratio of 
PLFA concentration in fungi to PLFA concentration in bacteria was sensitive to extreme 
precipitation events. The abundance and composition of soil bacteria were sensitive to extreme 
precipitation events. Alternatively, the ratio of PLFA concentration in gram-negative bacterial 
biomass to PLFA concentration in gram-positive bacterial biomass was used as a stress indicator 
(Klamer and Baath, 1998), as gram-negative bacteria have a thicker cell wall than gram-positive 
bacteria. Therefore, gram-positive bacteria are more resistant to drought than gram-negative 
bacteria. However, the ratio of PLFA concentration in gram-negative bacterial biomass to PLFA 
concentration in gram-positive bacterial biomass was not influenced by extreme precipitation 
events in our research. Compared with the moderate precipitation events, the extreme precipitation 
events did not significantly decrease SWC, although SWC had been shown to have contributed to 
the soil microbial community structure (Tiemann and Billings, 2011; Zeglin et al., 2013). 

At the phylum level, Actinobacteria dominated changes in soil bacterial composition under 
extreme and moderate precipitation events (Fig. 7c). The mean relative abundance of 
Actinobacteria was 43.8% higher under the P6 treatment than under the P24 treatment. The 
subdominant groups of dominant phylum were Proteobacteria and Acidobacteria. As previous 
studies showed that, under the simulating extreme drought treatments, the relative abundance of 
Actinobacteria (an oligotrophic bacteria phylum) increased, while that of Acidobacteria (a 
copiotrophic bacteria phylum) decreased (Fierer et al., 2007; Barnard et al., 2013). The different 
responses of soil microbial communities were due to the changed soil nutrient pool (Barnard et al., 
2013). The gram-positive bacteria (such as Actinobacteria) with high guanine and cytosine 
content were able to acquire nutrition interest as adverse conditions emerged (Zvyagintsev et al., 
2007). In many studies, Actinobacteria was found to be more resistant to drought than 
Acidobacteria (Goodfellow and Williams, 1983; Zvyagintsev et al., 2007; Barnard et al., 2013). 
Similarly, the abundance of Acidobacteria decreased with increased amount of precipitation 
(Zhang et al., 2013). In phylum Firmicutes, the organisms with gram-positive cell walls and low 
guanine and cytosine content produced endospores that are resistant to desiccation and thereby 
they can survive in extreme conditions. This partly explained why Firmicutes was only found 
under the P6 treatment. 

In general, the extreme precipitation events increased soil bacterial abundance, even reaching 
up to 40%, compared with the moderate precipitation events. However, we found no significant 
changes in soil microbial community structure and biomass among the precipitation events. This 
may be more related to the adaption strategy of soil microbe to extreme climate, such as extreme 
precipitation. Different soil microbial phyla have different life strategies (Fierer et al., 2007; 
Barnard et al., 2013). In addition, soil moisture availability has long been a major constrained 
factor to ecosystem productivity in the grasslands of Inner Mongolia (Bai et al., 2004; Hao et al., 
2010). Thus, soil microbial community may have adapted to the changes of precipitation patterns 
and increased resistance to extreme precipitation events (Gong et al., 2004; Bouskill et al., 2013). 


5 Conclusions 


The present study has important implications for understanding and predicting the impacts of 
changed precipitation events (large in amount and low in frequency) on soil microbial community. 
Our results illustrated the importance of extreme precipitation frequency, not merely precipitation 
amount, on soil microbial community structure and total soil microbial biomass. We found that 
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extreme precipitation events increased soil bacterial abundance, compared with moderate 
precipitation events. However, soil microbial community structure, total soil microbial biomass 
and soil fungal abundance remained quite stable in the extreme precipitation events. Soil bacteria 
were sensitive to extreme precipitation events in terms of numbers and classifications. These 
results implied that soil microbial community in the semi-arid steppe of the Xilin River Basin 
responds differently to changes in precipitation patterns. Soil bacteria and fungi may have 
different strategies to the changed precipitation events. We recommend long-term experiments on 
the influences of extreme precipitation events on soil microbial community and ecosystem 
functions across different ecosystems in the future. 
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